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Astrophysical neutrinos are powerful tools to investigate the fundamental properties of particle 
physics through their flavor content. In this paper, we perform the first general new physics study 
on ultra high energy neutrino flavor content by introducing effective operators. We find that at the 
current limits on these operators, new physics terms cause maximal effects on the flavor content, 
however, the flavor content at Earth is confined to a region related to the assumed initial flavor 
content. Furthermore, we conclude that a precise measure of the flavor content at Earth will 
provide orders of magnitude improvement on new physics bounds. Finally, we discuss the current 
best fits of flavor content of the IceCube data and their interplay with new physics scenarios. 
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Introduction — The existence of extra-terrestrial ultra 
high energy neutrinos has been confirmed by the Ice- 
Cube neutrino observatory m, opening the possibility 
to study ultra high energy particle production mecha¬ 
nisms as well as new neutrino physics 13 ®- The nature 
of these neutrinos from 35 TeV to 2 PeV is still a puz¬ 
zle, at the moment there are many astrophysical and be¬ 
yond the standard model candidate sources EHH2 that 
may produce these neutrinos. Currently, there is no sta¬ 
tistically significant spatial correlation between observed 
neutrinos and potential sources mm- 

Even though the sources of these neutrinos remains un¬ 
known, it is still possible to find evidence of new physics. 
The vacuum neutrino propagation Hamiltonian is lin¬ 
early proportional to the neutrino square mass differ¬ 
ences and inversely proportional to the neutrino energy. 
For astrophysical ultra high energy neutrinos this oper¬ 
ator is suppressed allowing to look for extremely tiny 
new physics effects otherwise cannot be seen. In the 
standard oscillation scenario, for any given initial flavor 
composition, the final composition, after the propaga¬ 
tion, lies in a small region on the flavor triangle close to 
(4> e : (ji^ : <f> T ) = ( 1:1:1). The flavor content of the 
astrophysical neutrinos has been studied in mmmm- 
These analyses find flavor content is statistically consis¬ 
tent with the standard oscillations expectations. Future 
data will clarify the IceCube astrophysical event flavor 
composition. 

In this paper, we perform the first general new physics 
study of the astrophysical neutrino flavor content by in¬ 
troducing effective operators in the standard three neu¬ 
trino scenario with unitary evolution. This is the far most 
general approach to study new physics in astrophysical 
neutrino flavors, and this approach covers many exotic 
particle physics models. There are few cases we do not 


consider in this paper. First, the model we work are lim¬ 
ited within lepton number conservation, and we do not 
consider models such as the neutrino-antineutrino oscilla¬ 
tions Hang. Second, we do not consider neutrino decay 
model which violates unitary evolution and was discussed 
elsewhere [20]. Similarly, we also do not consider mod¬ 
els with sterile neutrino states m- The sterile neutrino 
mixing matrix elements are known to be miniscule com¬ 
paring with the active neutrino mixing elements [211425] , 
and the contribution to the transition probability due to 
the sterile neutrinos is suppressed by the sterile-active 
matrix element to the fourth power. 

Ultra High Energy Astrophysical Neutrino Oscillations 
- Neutrinos change lepton flavors as they propagate 
macroscopic distances. This is due to the fact that the 
neutrino propagation eigenstates are not the eigenstates 
of the charged current weak interaction. In presence of 
a dense medium the decoherent scattering interactions 
are important [26], but in this paper we assume vacuum 
propagation. 

In general the relation between the propagation eigen¬ 
states | vf), and the flavor eigenstates \v a ), is given by a 
unitary transformation V(E), 

Wa) = Y^ Vai ( E )Wi) ■ M 

i 

For astrophysical neutrinos, the propagation distance is 
much longer than the oscillation length, and in this limit 
the oscillation from flavor state \v a ) to a flavor state \va) 
can be averaged, 

Pu a ^ p (E) = ^ \V a i(E )\ 2 \Vpi(E )\ 2 , (2) 

i 

where the probability depends only on the mixing matrix 
elements \V a i(E)\, which is in general energy dependent. 
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Using the probability given in this equation and the flux 
at production we can calculate the neutrino flux at 
Earth, for a flavor /3. It is more convenient to 

define the energy averaged flavor composition as 

Y. P ^^E)r a {E)dE , (3) 

Jae „ 


where we assume E~ 2 power law for the production flux 
and A E =[10 TeV, 10 PeV]. Note, however, that our 
main results are largely insensitive to the spectral index. 
We also assume that all flavors have the same energy 
dependence at the source. 

In astrophysics charged pion decay from proton-proton 
collisions is one of the preferred neutrino production 
channels. In this scenario the initial flavor composition 
is (</> e : (j)^ : 4> t ) = (1 : 2 : 0). Other scenarios such as 
rapid muon energy loss produce (0:1:0), neutron de¬ 
cay dominated sources produce (1:0:0) are of interest, 
while compositions such as (0 : 0 : 1) are not expected in 
the standard particle astrophysics scenarios. In order to 
plot the flavor content in a flavor triangle we introduce 
the flavor fraction, a ® = (j>® / 

For the vacuum propagation, the Hamiltonian of the 
standard neutrino oscillation only depends on the neu¬ 
trino mass term, 


1 f 0 0 0 

H= —U 0 Aro^ 0 

y 0 0 Am|-| 


U t = —UM 2 tf , (4) 


where E is neutrino energy, A mjj = m 2 — m|, and U 
is the standard lepton mixing matrix U. Throughout 
this paper we assume the normal mass ordering. We 
also performed the same study by assuming the inverted 
mass ordering, however, differences are minor and mass 
ordering does not affect any of our main conclusions. 

The current measurements of the standard neutrino 
oscillation experiments allows us to determine the astro- 
physical neutrino flavor content at detection given an as¬ 
sumption of the neutrino production. In Fig. [l] we show 
allowed regions of the flavor content at Earth, where we 
use the standard mixing angles and their errors from the 
global fits [27j in order to produce probability density 
distributions for the flavor content. Since the CP-phase 
is not strongly constrained by neither terrestrial [25] [25] 
nor astrophysical m neutrinos, we assume a flat distri¬ 
bution from 0 to 27 t. Note that for simplicity we use 
the larger of the asymmetric errors and implement them 
as Gaussian. In the left plot, we assume four different 
production flavor composition hypotheses. We observe 
that all the allowed regions of astrophysical neutrino fla¬ 
vor content at Earth are close to (1 : 1 : 1), except when 
the initial flavor content is (1 : 0 : 0) [31]. In the right 
plot, we show the allowed region of the flavor content of 
the astrophysical neutrinos with all possible astrophys¬ 
ical production mechanisms, i.e., the production flavor 



FIG. 1: Allowed regions of the flavor content at Earth using 
the priors on the mixing angles and errors given from the cur¬ 
rent neutrino oscillation measurements. In the left plot, the 
different colors correspond to different assumptions on flavor 
content at the production. The color intensity is proportional 
to the probability density. In the right plot, we further sample 
the initial flavor content as (x : 1 — x : 0). 


composition is sampled with (x : 1 — x : 0) uniformly on 
x [32] , Therefore, this rather narrow band covers all pos¬ 
sible scenarios of the standard neutrino oscillations with 
the standard astrophysical neutrino production mecha¬ 
nisms. 

New Physics in Effective Hamiltonians — An effective 
way of introducing new physics in neutrino oscillations 
is by introducing new operators. The full Hamiltonian 
that incorporates the new physics operators, in the flavor 
basis, can be expressed as, 

H = 2^ M2 tf f +£ (£) U n O n Ul = V\E)AV{E) , 

where O n = diag(O n ,i, 0„ i2 , O n , 3 ) and A = 
diag(Ai, A 2 , A 3 ). O n and A n set the scale of the new 
physics and U n is the mixing matrix that describes the 
new physics flavor structure. In the effective theory ap¬ 
proach, lower order operators are more relevant, thus in 
this work we will only study the first terms in the expan¬ 
sion, namely n = 0 and n = 1. 

Although in this work we will study n = 0 and n = 1, 
results can be extended to higher orders. These new oper¬ 
ators can be interpreted in different new physics contexts. 
Some examples for n = 0 new physics are couplings be¬ 
tween neutrinos and space time torsion [331, CPT-odd 
Lorenz violation [34H37] . and non-standard neutrino in¬ 
teractions [38MT1 . As for n = 1 new physics operators, 
CPT-even Lorentz violation [42}[43| and equivalence prin¬ 
ciple violation [21S5] are possible examples. 

There are some constraints from neutrino oscillation 
experiments to these effective operators in the context 
of Lorentz and CPT violation [46]. The most stringent 
limits on certain parameters are obtained from Super- 
Kamiokande and IceCube atmospheric neutrino analy¬ 
ses '!?, 25]. In this context the CPT-odd and CPT- 
even Lorentz violation coefficients are constrained to be 
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FIG. 2: Allowed region using anarchic sampling on the mix¬ 
ing angles for the new physics operator when the mass term in 
the Hamiltonian is neglected. The different plots correspond 
to different assumption on flavor content at production. The 
color intensity is proportional to the probability predicted by 
anarchic sampling. 


FIG. 3: Allowed region using anarchic sampling on the mix¬ 
ing angles for the new physics n = 0 operators. The top plot 
corresponds to the current limits on n = 0 operator; the bot¬ 
tom left plot corresponds to O o = 3.6 x 10 -26 GeV, while the 
bottom right plot corresponds to Oo = 6.3 x 10 -28 GeV. 


~ 10~ 23 GeV and ~ 10 -27 depending on the flavor struc¬ 
ture U n . These constraints can be used to set the scales 
of n = 0 and n = 1 operators introduced in this paper. 
For example, we set Oo = 1 x 10 -23 GeV as a current 
limit of the n = 0 operator, and 0\ = 1 x 10 -23 GeV with 
Ai = 1 TeV as a current limit of n = 1 operators, where 
= 10 -27 . Through this paper we have assumed the 
scale of 0\ is of the order of O 0 without loss of generality. 

Anarchic Sampling prediction and IceCube Results — 
In order to predict the flavor composition at Earth in the 
presence of new physics, the values of the mixing matri¬ 
ces U n should be specified. In order to show a prediction 
with new physics operators, we have to account for all 
the free parameters in the mixing matrix; we use a ran¬ 
dom sampling scheme to construct the mixing matrix. 
A well established schema is the anarchic sampling m- 
13~2j . which samples a flat distribution given by the Haar 
measure, 

dU n = ds\ 2 A dc\ 3 A ds\ 3 A dS , (5) 

where, , Cij, and S correspond to sines, cosines, and 
phase for the new physics n-operator mixing angles. We 
omit the Majorana phases since they do not affect neu¬ 
trino oscillations. 

In Fig(2] we show the allowed regions using anarchic 
sampling in the case where H = ) U n O n U\. In this 


case, we neglect the mass term and we are considering 
that the Hamiltonian has only one operator, i.e., V = 
U n , and the result does not depend on n. Each plot 
in this figure correspond to a different production flavor 
composition. We show the pion decay production (1:2: 
0) [yellow], beta decay (1:0:0) [green], muon cooling 
(0:1:0) [red] and for completeness we show the exotic 
v T dominant model (0:0:1) [blue]. The color density 
in these plots is a representation of the probability given 
by the anarchic sampling. 

In Figj3]we show the case where we have a mass term 
and the n = 0 operators. In the top plot, we set Oo = 
l.OxlO -23 GeV, corresponding to the order of the current 
best limit on this operator. On the bottom left plot we 
set Oq = 3.6 x 10 -26 GeV and the bottom right plot 
we set Oo = 6.3 x 10 -28 GeV. These values are chosen 
because they have the same magnitude as the mass term 
with neutrino energy of E v = 35 TeV and E v = 2 PeV 
respectively. In this plot, the colors represent different 
assumptions in the production flavor content, and the 
color intensity is the probability given by the anarchic 
sampling as in Fig [2] 

In Fig|4] we show the case for the n = 1 operators. 
The color notations and their intensities have equivalent 
meaning as Fig(3] As before, the top plot we set the new 
physics operator to the current best limit ^ ~ 10 -27 . 
This is achieved by choosing 0\ = Oo = 1.0 x 10 -23 GeV 











4 



FIG. 4: Allowed region using anarchic sampling on the mix¬ 
ing angles for the new physics n = 1 operators. The top 
plot corresponds to the current limits on n = 1 operator; the 
bottom left plot corresponds to Oi = 3.6 x 1CP 26 GeV and 
Ai = 35 TeV = 1.0 x 10~ 3 °), while the bottom right 
plot corresponds to Oi = 6.3 x 10 -28 GeV and Ai = 2 PeV 
(gi = 3.2 x lO" 34 ). 


and Ai = 1 TeV. In the bottom left plot, 0\ = 3.6 x 
10~ 26 GeV and Ar = 35 TeV are used, and in the bot¬ 
tom right plot the parameters are 0\ = 6.3 x 10~ 28 GeV 
and Ai = 2 PeV. These choices make new physics to 
be the same magnitude as the mass term with neutrino 
energy of E v = 35 TeV and E v = 2 PeV, respectively. 
In other words, these choices explore new physics down 
to ^ = 1.0 x 1CT 30 and ^ = 3.2 x 10" 34 . This can be 
compared, for example, to the aforementioned best limits 
of Lorentz and CPT violation in neutrino sector 03 mi. 
The potential limits from astrophysical neutrino flavor 
content can be well beyond what terrestrial neutrino ex¬ 
periments can achieve. 

From Fig(3]and Fig|4]we observe that the allowed re¬ 
gions in the flavor triangle change in a similar way as a 
function of the energy scale. This is true for any higher 
operators, because what matters is the scale where they 
dominates over standard neutrino mass term, and these 
two operators are sufficient to predict behaviors of any 
higher order operators. Comparing Fig(3]and Fig|4]with 
respect to Fig(2]where the allowed regions are more sym¬ 
metric, there is a preferred region along the vacuum oscil¬ 
lation triangle shown in FigjT] It is interesting to notice 
that due to the unitary evolution and the fact that the 
oscillations are averaged, for a given production flavor 


content, only a subset of the flavor triangle is accessi¬ 
ble. The pion decay production mechanism (1 : 2 : 0) is 
one of the most natural astrophysical scenarios for high 
energy neutrino production. From Fig(3] and Fig|4] the 
allowed region for this case is the smallest, which means 
that if future measurements exclude this region, the pion 
production dominant mechanism is excluded regardless 
of the presence of new oscillation physics. 

In the analyses of the IceCube high energy neutrino 
events, different results have been shown. The first re¬ 
sult 03 ], using the IceCube result 0 ], showed a best fit 
at (1 : 0 : 0) disfavoring (1:1: 1) at 92% C.L. Later, the 
same authors did an improved analysis 01 including en¬ 
ergy dependence and extra systematic errors, finding that 
the best fit may move considerably depending on the fea¬ 
tures of the energy spectrum such as including an energy 
cutoff or not. The IceCube collaboration later published 
an analysis of the flavor ratio above 30 TeV m finding a 
best fit at (0 : | : |), as well as excluding (1:0:0) and 
(0 : 1 : 0) at more than 90% C.L. This IceCube result 
shows a best fit dominated by the v T component, which 
can be explained by the correlation between the energy 
cutoff and the Glashow resonance, as noted by [14]. In 
obtaining this best fit, the IceCube collaboration has as¬ 
sumed an equal amount of neutrinos and antineutrinos, 
which best corresponds to a proton-proton source. On 
the other hand, if the neutrino source is proton-photon 
dominated then the neutrino-antineutrino ratio weaken 
making the previous conclusion. It is interesting to notice 
that if this IceCube best fit does not change considerably 
after adding more data, the production mechanism has to 
include a v T component. This is because the new physics 
in the propagation can not give the best fit value for any 
plausible astrophysical scenarios. This implies not only 
new physics in the neutrino oscillations, but also new 
physics in the production mechanism. 

Conclusions — We performed the first new physics 
study on the astrophysical neutrino flavor content using 
effective operators in the standard three neutrino sce¬ 
nario. These operators can represent a variety of models 
such as Lorentz and CPT violation, violation of equiva¬ 
lent principle, cosmic torsion, non-standard interactions, 
etc, making this work to be the most general study of 
new physics in astrophysical neutrino flavor content to 
date. 

We found that large effects in the flavor content at 
Earth are still allowed with given terrestrial bounds on 
new physics in the neutrino sector. This implies that an 
accurate measurement of the flavor content will provide 
stronger bounds on new physics. Furthermore, there are 
regions on the flavor triangle that cannot be accessed 
even in the presence of new physics in the neutrino os¬ 
cillations for any of the plausible astrophysical mecha¬ 
nisms. Interestingly, that the most natural astrophysical 
mechanism, pion decay, has the smallest region in the 
flavor triangle even when new physics is considered. The 
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real astrophysical neutrino production mechanism in na¬ 
ture may be the combination of channels, but our results 
hold for such a case. Therefore, a higher statistics mea¬ 
surement by future neutrino telescopes, such as IceCube- 
Gen2 [M] , could reveal not only the initial neutrino flavor 
ratios, but also the presence of new physics in neutrinos. 
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